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Abstract. In this paper, we first obtain an algebraic formula for the moments of 
a centered Wishart matrix, and apply it to obtain new convergence results in the 
large dimension limit when both parameters of the distribution tend to infinity at 
different speeds. 

We use this result to investigate APPT (absolute positive partial transpose) 
quantum states. We show that the threshold for a bipartite random induced 
state on C'' = (g) C^^, obtained by partial tracing a random pure state 
on (E) C''*, being APPT occurs if the environmental dimension s is of order 
So = min((ii, ^2)'^ max((ii, (i2). That is, when s ^ Csq, such a random induced 
state is APPT with large probability, while such a random states is not APPT 
with large probability when s ^ csq. Besides, we compute effectively C and c and 
show that it is possible to replace them by the same sharp transition constant when 
min((ii, (i2)^ < d. 



1. Introduction 

Geometry of quantum states strives to understand the geometric properties of sub- 
sets of quantum states, and has attracted considerable attention, especially in the 
case of the large dimension [2, 3, 15, 35, 37, 38, 41, 42]. The high dimensional set- 
ting is common (and of particular interest) in Quantum Information Theory whose 
building blocks, quantum states, often are objects with huge dimension (for instance, 
the set of quantum states on the system (C^)®^ has dimension 43046720). This high 
dimensional setting indicates the importance of random constructions which now is 
a main tool in understanding the typical behavior of random induced states. To gen- 
erate random induced states, one often relies on random matrices. The connections 
between Random Matrices and Quantum Information Theory were pushed forward 
by Hayden, Leung, and Winter in their studies of aspects of generic entanglement 
[19]. Together with tools from Geometric Functional Analysis and Convex Geometric 
Analysis, random matrices and random constructions have led to many important 
(and even unexpected) results, such as Hastings's disproof of the famous additivity 
conjecture for the classical capacity of quantum channels [18]. Recent contributions 
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include the studies of the generic properties for entanglement vs. separability, and 
PPT (positive partial transpose) vs non-PPT [1, 4, 5, 10]. 

Detecting quantum entanglement, a phenomenon first discovered in [13] and now 
being the key ingredient of quantum algorithms (see [27, 32]), is one of the funda- 
mental problems in Quantum Information Theory. Among those necessary and/or 
sufficient conditions for separability and entanglement, the Peres- Horodecki PPT 
criterion [22, 30] is the simplest but the most powerful one. The Peres-Horodecki 
PPT criterion is a necessary condition and is sufficient only for the systems ® 
and ® C'^ [34, 40] . From the computational complexity point of view, separability 
and PPT are quite different: determining separability is an NP-hard problem [16], 
but determining PPT is easy since it only requires to verify the eigenvalues of the 
partial transpose of given states being positive. Note that both the separability and 
PPT are encoded in the spectral properties of quantum states in a complicate way. 
Necessary and/or sufficient conditions on determining separability and PPT by just 
the information of eigenvalues (referred to as the absolute separability and absolute 
PPT in literature) could be very useful in Quantum Information Theory as it can 
help to reduce the cost of storage spaces and (processing) time. For absolute separa- 
bility, less results are known; however, necessary and sufficient conditions for APPT 
have been found by Hildebrand [20]. Understanding when a random induced states 
is APPT is the main motivation of this work. 

To that end, we first prove that a properly centered d x d Wishart matrix of pa- 
rameter s has its expected normalized moments that can be written as a polynomial 
in the variables d and d/s. The coefficients of this polynomial have a simple combina- 
torial interpretation, and some families of coefficients are known. This algebraic fact 
has important consequences in the two parameter asymptotic study of the Wishart 
matrix. Indeed, it allows in particular to capture the precise nature of the behavior 
of the Wishart matrix in the case where d ^ oo, and in particular in the case where 
s/d— j-ootoo. We summarize our first main result as follows: 

Theorem A. Let Wd be a d x d Wishart matrix of parameter s and let Zd = 
\fds (-^ ^ ^) be its centered and renormalized version. The moments of Zd are 



Moreover, almost surely as d —)■ oo and s/d ^ oo, the extremal eigenvalues of Zd 
converge to ±2. 

Both results (for moments and for extremal eigenvalues) are of separate interest in 
random matrix theory, where the single scaling d — oo, s/d — )■ c > has received a 
lot of attention[14, 23, 25]. The above result is then applied to estimate the threshold 
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for a random induced state being APPT vs. non-APPT. We have the following result 
(we put p = min((ii, 6/2)): 

Theorem B. There are effectively computable absolute constants c,C > 0, such 
that, if p is a bipartite random induced state on = C^^ ® C'^'^ , obtained by partial 
tracing over C a random pure state on (g> C*, then for d = did2 large enough, 
one has: 

(i) The quantum state p is APPT with very large probability when s ^ {A + e)p'^d. 

(ii) The quantum state p is not APPT with very large probability when s ^ cp'^d. 
If p^ <^ d, one can take c = 4 — e; when p^ is of the order of d, c can be 
computed as described in the proof of Theorem 4-i- 

The letters C, c, cq,... denote absolute numerical constants (independent of any- 
thing) whose value may change from place to place. When A, B are quantities de- 
pending on the dimension (and perhaps some other parameters), the notation A < B 
means that there exists an absolute constant C > such that the inequality A ^CB 
holds in every dimension. Similarly A ^ 5 means both A < B and B < A. As usual, 
A B means that A/i? — > 1 as the dimension (or some other relevant parameter) 
tends to 00, while A = o{B) means that A/B — j- 0. For a d x d complex matrix 
A G Mrf(C) we denote by tr(y4) its non-normalized trace. In this paper, whenever 
we deal with a tensor product structure = C^^ (g) C^^, we put p = min((ii, ^2)- 

The article is organized as follows. Section 2 is of random matrix theoretic flavor. 
Using mostly combinatorics and also a little bit of analysis and elementary probability 
we obtain new formulas for the moments of centered Wishart matrices and estimates 
on their extremal eigenvalues. Section 3 gathers some properties about the the APPT 
property and section 4 provides the bounds of the threshold for APPT. 

2. Combinatorics of centered Wishart matrices 

In this section, we prove two results. Theorems 2.2 and 2.7 about centered and 
renormalized Wishart matrix. These results are interesting for random matrix the- 
orists and can be considered independently from the rest of the paper. 

2.1. Preliminaries and notation. We start by introducing some notation from 
combinatorics. For an integer p, we denote \p] = {1,2, ... ,p}, and [0] = 0. For 
a subset / C [p], let Sj be the set of permutations which act on /. We shall take 
the convention that 5© = {0}. The set of permutations without fixed points will de 
denoted by 

S'^ = {a e Si \yi e I, a{i) ^ i}. 

The length \a\ of a permutation a G 5/ is the minimal number of transpositions 
needed to decompose a. We put |0| =0. The notation | ■ | is polymorphic, since it is 
used to denote both the cardinality of sets and the length of permutations. We shall 
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also use the notation jj^a for the number of cycles of a; the following relation holds 
#a + |a| = |/|. 

For a nonempty subset J, we denote by 7/ the full cycle in 5/, with elements of 
/ ordered increasingly: 77 = (zi Z2 ■ ■ ■ ^fc), where I = {ii < 12 < ■ ■ ■ < ik}- Abusing 
notation, we define I70I = — 1. In this way, the geodesic inequality 

|a| + |a~"^7/| ^ I7/I = |/| — 1 

holds for all I and all a G 5/. We define the genus of a permutation a as half of the 
amount by which the above inequality fails to be an equality 

, , lal + |a~"^7/| -1/1 + 1 
^/(«) = ^ • 

It is a standard fact in combinatorics that the genus gi{a) is a nonnegative integer. 
For a G iS/, define a G iSj to be a without its fixed points; in other words, J = {z G 
/ I a{i) 7^ 1} and, for j G J, we have a{j) = a{j). It is clear that \a\ = \a\. Moreover, 
by the following lemma, erasing fixed points leaves the genus of the permutation 
unchanged. 

Lemma 2.1. Let a E Sj be the permutation a E Sj with its fixed points removed. 
Then, gj{a) = gi{a). 

Proof. Going in the opposite direction and proceeding by induction, it suffices to 
show that whenever we add a fixed point i to a permutation a G iSj, its genus 
remains unchanged. Let us denote by ji and j2 the neighboring points in J between 
which i is inserted: ji < i < j2. Also, we note by ji — 1 the predecessor of ji in 
J and by j2 + 1 the successor of j2 in J- Using the number of cycles notation, one 
needs to show that #(«~^7j) = #(a~^7/). 

Given two permutations a, r G Sp, recall the following combinatorial interpretation 
of the number of cycles of Define a multigraph Go-,t = {V, E) with 2p vertices 

V = {1, . . . ,p,l' , . . . ,p'} and edges 

E = {ik, aik)') ,ke[p]}U{ik, T{k)') , k G [p]}. 

Then the number of cycles 7^((T^^r) equals the number of connected components of 

Going back to our setting, it is clear that the graphs Ga,7j and Go,-y^ have the 
same number of connected components, since adding the extra vertices i,i' does not 
alter the edge structure of Ga^^j] for a sketch of the argument, see Figure 1. □ 

2.2. A moment formula for the centered Wishart matrix. Let G G Mdxs(C) 
be a Ginibre random matrix (i.e. {Gjj} are i.i.d. standard complex Gaussian random 
variables) and W = = GG* be the corresponding Wishart matrix of parameters 
{d,s), where G* denotes the Hermitian adjoint of G. Here we make an abuse of 
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Figure 1 . Adding a fixed point i to a permutation a does not increase 
the number of cycles of a'^jj. Edges corresponding to a and a are 
represented by dasfied red lines and edges corresponding to full cycles 
are black and solid. 

notation and keep track only of the parameter d, considering implicitly that s will 
be a function of d. It is easy to see that KWij = s6ij and EtrlV = ds. 

The main theorem of this section characterizes the fluctuations of W around its 
mean. For this purpose, we introduce the following dxd matrix, a properly rescaled, 
centered Wishart matrix of parameters {d, s): 



where 1^ refers to the dxd identity matrix. We show that the following theorem 
holds true. 

Theorem 2.2. The moments of Zd are given by 



Note that despite the simplicity of this combinatorial result, it seems to be new. 
Before we prove this result, we would like to describe three corollaries, obtained by 
letting one or both parameters d and s go to infinity. 

Corollary 2.3. If d is fixed and s — t- 00, then 




(1) 





if p is odd, 
if p is even, 
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where e{p/2,g) is the number of products ofp/2 disjoint transpositions in Sp of genus 
g. Alternatively, for even p, e{p/2,g) is known to count the number of gluings of a 
p-gon into a surface of genus g. 

Proof. When s— 7-00,^/5— ^-O the only terms in equation (1) which survive are those 
for which |a| = p/2. It follows that a must be in this case a product of p/2 disjoint 
transpositions (for even p). Reordering the sum by genera gives the statement (see 



Note that the sequence e{p,g) appears in the Online Encyclopedia of Integer Se- 
quences [29] as A035309. 

For the forthcoming corollary, we need to recall that the set NC{p) is the collec- 
tion of partitions of [p] that have no crossings with respect to the canonical order. 
Moreover, we introduce the subset 



Corollary 2.4. // both rf, s — 00 such that s/d — c for some constant c > 0, we 
obtain 



where ^ir denotes the number of blocks of the partition vr. In particular, the random 
matrix Zd converges in moments to a centered Marchenko-Pastur distribution of 
parameter c (rescaled by \fc). 

Proof. In this asymptotic regime, the surviving terms in equation (1) are those for 
which gia) = 0. The formula in the statement follows from a well known result 
of Biane [8] saying that the permutations in Sp of genus are in one to one corre- 
spondence with non-crossing partitions vr G NC{p). The second part follows from a 
centered version of the free Poisson limit theorem [26, Theorem 12.11]. For c = 1, 
the rescaled quantities appearing in the statement are the Riordan numbers Rp [29, 
sequence A005043 ] such that Rp = \NC°{p)\; see [28] for the connection between 
Riordan numbers and centered free Poisson random variables. □ 

Finally, we have the following general asymptotics 

Corollary 2.5. If d —> 00 and s/d ^ 00 (in other words 1 d <^ s), then 



where Cat„ is the n-th Catalan number. In particular, the random matrix con- 
verges in moments to a standard semicircular distribution. 



[43]). 



□ 



NC°{p) = {tt G NC{p) I vr has no singletons}. 





if p is odd, 
if p is even., 
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Proof. This follows directly from the fact that the Catalan numbers count the number 
of non-crossing pairings of [2p\: e{2p,0) = Catp. □ 

Note that the above result also follows from the more general result by [7], where 
the almost sure convergence is also obtained. A proof of the almost sure convergence 
could also be obtained in the combinatorial spirit of this paper, for instance along 
the lines of [12]. 

We would like to explain briefly why this result is not surprising from a heuristic 
point of view. Indeed, the distribution for Marcenko-Pastur distribution, as given in 
Proposition 4.3, is 



^/4c - (x - 1 - c)2^ 
27rx 

so as c — oo, the distribution of Zd should approach 



, J (2 + - - 2 + ^/dfs) 

s/d— ; Ir 



\/s/d— ; If „, /TTn, /TrX^)dx 

which should tend to the semi circle distribution. Our corollary therefore implies 
that we let s / d go to infinity separately or together with any correlation we like. 
Let us now prove the combinatorial result. 

Proof of Theorem 2.2. The starting point is a formula for the moments of W , ob- 
tained via the Wick calculus (for an intuitive graphical approach to this problem, 
see [12]): 

adSp adSp 

where 7 G iSp is the forward cycle 7 = (12 ■ ■ ■ p). By applying the binomial formula 
for the commuting matrices W and 1^, we get 

ds d 

d-^+p/^sP'^Y^-lf'\dsf\d-\''\mi{W\'\) 

Y Y (-l)P^I^I(i^l-f/2+|/h|a-l7lsP/2-|a| 

\a\-p/2 



nip := E— tr 
d 



J2 ^(-l)^'H/|rf-2,,(.) 



/C[p] aeSj 
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To conclude, we need to show that the terms in the sum above cancel out, except for 
the ones with I = [p] and a G S[p] without fixed points. For a permutation a G Sj, 
denote by [a] the set of permutations which extend a by adding fixed points: 

[a] = {a e Si \ J C I , a{j) = a{j) Vj G J and a{i) = i Vi G / \ J}. 

Regrouping terms in the sum and using the fact that for a G [a], \a\ = \a\ and 
gi{a) = gj{a), we can write 

\a\-p/2 

rrir. 



565° ^ ^ «G5/n[Q 

where the first sum in the equation above is indexed by permutations a without 
fixed points. Given such a permutation a G iSj, for every larger set / D J there is a 
unique way of extending a to a E Si. Hence, the second sum in the above equation 
is given by 

which can be understood as a Mobius inversion formula in the poset of subsets of [p] . 
In conclusion, only the permutations a G iS° give non-zero contribution. To finalize 
the proof, note that such a permutation has at most p/2 cycles and thus at least 
length p/2. □ 

2.3. Almost sure convergence of extremal eigenvalues of the Wishart ma- 
trices in the regime d <^ s. We start with the following lemma: 

Lemma 2.6. Let us assume that d/s c E (0, oo) as d ^ oo. For any e > there 
exists a constant C > such that for any p ^ y/d, we have 

EtrZ^ ^ C{2 + ^ + ey. 

Proof. Under the hypotheses of the lemma, the eigenvalues counting measure of Zd 
converges almost surely to a probability measure whose support is [—2 + ^/c, 2 + ^/c\. 

Let us first evaluate PdlZ^H ^ 2 + a/c), where \\Zd\\ is the operator norm of Z^. 
According to the union bound, this is bounded above by 

P(Ai(Z,) ^ 2 + v^) + nUZd) ^ -2 - v^), 

where Xi{Zd) and \d{Zd) are respectively the largest and smallest eigenvalues of Zd. 

According to Theorem 3 of Soshnikov in [33], there exists a constant C > such 
that for any p ^ Vrf, we have 



\2p 
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By Jensen inequality this implies under the same assumption on p that 

Besides, it follows from Equation (15) in [24], that the probability of having eigen- 
values less than t < {\/c— 1)^ is less than exp(— (i(yf(t)). For our purposes it is enough 
to know that g{t) > as long as < t < {y/c — 1)^. 

We can conclude the proof of the lemma from the two above observations via the 
inequalities 

and the formula 



/•oo 

E(X) = / P(X ^ t)dt. 
Jo 



□ 



Now we would like to let s/d — )■ oo simultaneously (but independently) with 
(i — 7- oo. This is the setting of our next result, namely: 

Theorem 2.7. Almost surely, when 1 <^ d <^ s, the extremal eigenvalues of Zd 
converge to ±2. 

Proof. Let e > and do large enough such that for d ^ do, \/d/s ^ e/2. The 
moment inequality from Lemma 2.6 together with the fact that our moment formula 
of Theorem 2.2 involves only positive terms implies that there exists a constant 
C > 0, such that for all p ^ y/d, 

E(trZ^) ^ C(2 + ef. 

Letting d — )■ oo, by the Borel-Cantelli Lemma and the same Jensen inequality as 
in the proof of Lemma 2.6, we obtain that limsup Ai ^ 2 + e almost surely. Since 
this holds true for all e > and since liminf Ai ^ 2 by Corollary 2.5, we get the 
desired result. □ 

An interesting aspect of the above proof is that it relies on moment techniques, 
and therefore makes use of Theorem 2.2. Here the interest of the moment method 
and combinatorics is that they explain why one can make rigorous a change of limit 
between d ^ oo and s/d ^ oo regarding the almost sure convergence of the largest 
eigenvalue. 

We could have obtained directly this result with complex analysis results (see 
e.g. [9]) and one could probably have obtained more refined estimates (e.g. large 
deviation bounds, universality results, etc); however, this was not in the spirit of our 
combinatorial approach and we leave it for future investigation. 
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Using similar techniques, one could extend the above results to show the following 
quantitative bound: for all e > 0, with exponential small probability in — )■ oo, the 
spectrum of a random density matrix p from the induced ensemble of parameters 
{d, s) is contained in the interval 

'1 _ 2(1 +£) 1 ^ 2(l + £) ' 

_d y/ds ' d y/ds 

3. Existence of a threshold for AWT 

3.1. Quantum states and Absolute PPT. We now introduce some necessary 
notation and concepts related to quantum information theory; readers are referred 
to [11] and [27] for more details. 

Consider a (complex) Hilbert space = C^^®- ■ ■®Q'^'' with (complex) dimension 
d = di ■ ■ ■ dk, where c/j ^ 2 for alH = 1, ■ ■ ■ , fc. The set of states on (denoted by 
V = ©(C^)) can be identified with the set of d x d density matrices. That is, 

V = {pe (C) I p ^ and trp = 1}. 

Clearly the real dimension of V is N = d"^ — 1. The set V is contained in the afiine 
hyperplane 

Hi = {Ae Mf (C) I tiA = 1}, 

endowed with the Hilbert-Schmidt inner product {A,B) = tT{A*B). 

Partial tracing states on Cg> C'^ over gives reduced density matrices of size 
d X d. Any state p on ® C'^ may be written as 

d s 

where Eiaji3 are the matrix units associated to orthonormal bases {ei}f^i and 

of and respectively. The partial trace of p over C^, denoted by a = trcs(p) 

may be formulated as 

s 

aij = ^ pifjjfj for i,j = 1, - ■■ ,d. 

13=1 

The induced measure on V^C^) by partial tracing over is an important proba- 
bility distribution and can be described as follows. Let be a random pure state 
on ® C'^, where ip is a. random unit vector uniformly distributed on the sphere 
in (8> C*. Then the random induced state p = trcd'?/') ("^M) G V^C^) follows the 
distribution pd,s- Equivalently, one can find a. d x s matrix M distributed uniformly 
on the sphere of d x s matrices, such that p = MM*. The distribution pd,s plays 
central roles in this section. When s = d, one gets pd,d, the normalized Lebesgue 
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(i.e., the Hilbert-Schmidt) measure on ©(C^). Hence, a random state distributed 
according to fid^d is uniformly distributed on V^C^). When s ^ d, the probabihty 
measure fid,.s has a simple form [44] 

(2) P^{p) = ^{detpr^ 

dvol Zd,s 

where Zd^s is a normalization factor. Note that formula (2) allows to define the 
measure pd,s (in particular) for every real s ^ d, while the partial trace construction 
makes sense only for integer values of s. 

Hereafter, we will focus on the bipartite system, i.e, the case C"' = ® C'^^ with 
d = did2- The partial transpose operator (denoted by is a linear operation that 
consists in taking the transpose in one leg and doing nothing on the other leg, i.e., 
T2{ti®T2) = Ti ®r^, where T is the normal transpose operator. The set of quantum 
states with positive partial transpose is denoted by WT, i.e., p G WT if and only 
if T2(p) ^ 0. Geometrically, WT = V (1 T2(T>), and WT is a convex body with 
constant height [36]. Peres- Horodecki PPT criterion states that S C WT [22, 30], 
and S = WT only if d ^ 6 [34, 40]. Here, the = rf^ - 1 dimensional set 5 C I' 
is the set of separable quantum states on [39] defined as 

5 = 5(^^1 ^C^^) :=conv{pi®p2, PiEViC"'), p2eV{C''')}. 

(Similarly, one can define 5(C'^i ® ■ ■ ■ ® C^").) The set S := V \ S is the set of 
entangled quantum states, which play a crucial role in quantum information and 
quantum computing. 

A quantum state p G ©(C^^ (g) C^^) jg absolutely PPT{ot APPT) if for any unitary 
matrix U G U{d), UpU* G WT. The set of all states being APPT, denoted as 
AWT, is 

(3) AWT = Pi U{WT)U* c WT. 

ueuid) 

Clearly, AWT is a convex body, a convex compact set with non-empty interior. This 
follows easily from (3) and the following result in [17]: + (1 — e)^ C 5 C WT 
for some e < ^zj- 

In applications, one often requires the convex bodies to be origin-symmetric. 
In this section, we will mainly work on the symmetric convex body AWTsym = 
-AWTo n AWTo where AWTo = AWT - l/d. Such a symmetrization will 
not change many geometric parameters of interest (such as, the volume radius and 
mean width) substantially, due to the famous Rogers- Shephard inequality [31]. Both 
AWTo and AWTsym sit in the linear hyperplane 

Ho = {Ae Mf (C) I tiA = 0}. 
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Recall that the gauge associated to a convex body K C R is the function || ■ Un- 
defined by 

\\x\\k = inf{t ^ : X e tK}, Vx G R^. 
Note that \\x\\k = || — x\\k for origin-symmetric convex bodies K. The outradius of 
K is the smallest R > such that K is contained in a ball of radius R. Similarly, 
the inradius of a convex body K is the largest radius r of a Euclidean ball contained 
in K. For origin-symmetric convex bodies, r and R can be defined as the "best" 
constants such that R~^ \ ■ | ^ || ■ ||k ^ '""^1 ■ | (where | ■ | denotes the Euclidean norm.) 
Let X be a standard Gaussian vector in R^, i.e., a random vector with independent 
A^(0, 1) coordinates in any orthonormal basis. 

Proposition 3.1. In the notation of the present section, we have 

(4) . I ^ II . m^^^^^^ ^d\-\ ^ d-'\ ■ K II ■ |Uppr°,„, ^ d-'/^\ ■ I, 

(5) E llxm^pro ^ E llxm^prsym ^ 2E Hxm^pro- 

Proof. Any matrix A G AWTq satisfies A ^ —I/d. This implies that any A G 
AVVTsym satisfies —I/d ^ A ^ I/d, or \\A\\oo ^ l/d, and therefore the outradius 
of AVVTsym is bounded by l/Vd. The inradius of AWT equal to {d{d - 1))-^/^ 
follows directly from Corollary 3 in [17]. This proves (4). 
Note that the distribution of X is symmetric, and 

||Ar||^Ppro ^ W^WAVVTsym 

= max(||Xmpp7-o5 II ~ ^WavvTo) ^ II-^IUppTo + II ~ ^WavvTo- 
Then (5) follows after taking expectation. □ 

It was pointed out that the set AWT varies if the decomposition of d changes: 
if min{d'i,d2) ^ min((ii, ^2), then any APPT quantum states on C^i must be 

APPT on (g) C^^ [20]. Consequently, the largest APPT set is obtained when 
di = 2 and d2 = d/2, and the smallest APPT set is obtained when di = d2 = \/d- 

Let p = min((ii, ^2)- Denote by 5+ = {{k,l) : 1 ^ k ^ I ^ p} and S'_ = {{k,l) : 
1 ^ k < I ^ p}. Note that the cardinahties of the sets 5*+ and S_ are p+ = 1)/2 
and p_ = p{p — l)/2 respectively. Let 

cr+ : 5+ — 7- {1, ■ ■ ■ and a"_ : 5- — t- {1, ■ ■ ■ ,P-}, 

be two orderings (i.e. bijective maps) on S+ and S'_ respectively. Thus, cr+(/c, /) ^ p+ 
and (J_{k, I) ^ p_ for all pairs (fc, /). 

For A = (Ai, ■ ■ ■ , \d) with ordering Ai ^ A2 ^ ■ ■ ■ ^ A^ ^ and Ylii=i -^j — 1' 
defines the matrix A(A; 0"+, cr__) as 
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Define the px p symmetric matrix 6(A; 0"+, cr_) to be the sum of A(A; cr+, (T_) and 
its transpose A-^(A; (t+, cr_). Note that 6 (A; (t+, cr_) has the following form: 

(6) ew=l 

' \ \b{k,i) - Kik,i), b{k,l) > c{k,l),k> I, 

where a{k), b{k,l) and c{k,l) are some integer- valued functions with values smaller 
than or equal to d. Thus Qk,i ^ for all k ^ I. 

The following theorem is a necessary and sufficient condition for p G AWT (see 
Theorem 111.9 or Lemma 111.10 in [20] ). 

Theorem 3.2. Let p he a quantum state on = C^^ CSC^^ with d = did2. Suppose 
that p has eigenvalues Ai ^ A2 ^ ■ ■ ■ ^ A^ ^ 0. Then p is AWT if and only if for 
all pair of ordering ((T+, cr_), 0(A; cr+, a ) is positive semi- definite. 

3.2. Threshold for AWT is ~ w{AWTlf. Hereafter, w{AWT°o) denotes the 
mean width of AWTq, the polar body of AWTq. For general convex body K 
with the origin in its interior, its polar body (denoted by K°) can be defined as 

K° = {ye'R^ : {x,y) ^ 1 Vx G K}, 

where (■, ■) denotes the usual inner product and induces the Euclidean norm | ■ |. The 
mean width of K, w{K), is defined as 



(7) w{K) := / hK{u)da{u)= / ||n||i^o(io-(n), 

where da{u) is the normalized spherical measure on the sphere S^~^, and hxiu) = 
max^.gi^(x, m) = ||m||/<o for any u G S^^^. A more convenient quantity to calculate 
is the Gaussian mean width of K 

(8) wg{K) =E\\X\\k- =E sup {X,x), 

where X is a standard Gaussian vector in R^. By passing to polar coordinates, one 
can easily check that for every convex body K C 

(9) wciK) =^Nw{K), 
where 

(10) y]v^«.„^V7v. 

is a constant depending only on N. We set SQ^di, ^2) to be 



, E\\X\Uvvro Y [ wg{AWTIj ^ ( a-v-vt^ 
So = So{di, 4) = ( — 1 = I — ) ~ w{AWTq 



•o\2 
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By inequality (5), one has Sq ~ w{AVVTl)^ - ^{AVVTlyJ^ . The following 
theorem states that the threshold for the set AWT is of order of w{AVVTqY . 

Theorem 3.3. There are effectively computable absolute constants C,c > 0, such 
that, if p is a random induced state on C^^ ^ C^^ distributed according to the measure 
fid,s, then 

(i) P(p e AWT) ^ Cexp(-cs) for s ^ csq; 

(ii) P(p ^ AWT) ^ Cexp(-cso) for s ^ Cso- 

We first point out that the threshold value for the set WT occurs only at those 
s ^ 2d. For the balanced bipartite case (i.e. di = 0^2) it follows from Theorem 4 in 
[1], while for the unbalanced bipartite case (i.e., di ^ 0^2) it follows from [10]. As 
AWT C WT, the threshold for AWT must be larger and thus we also have 
s ^ 2d. 

The following lemma is our main tool to prove that the threshold for the set AWT 
can be taken as w{AWTq)^. This lemma aims to approximate p — I/d hj -^^^X 
where X is a standard Gaussian vector in the space Hq of traceless Hermitian d x d 
matrices. We refer readers to its detailed proof in [4] . 

Lemma 3.4. For every convex body K C Hq containing in its interior, and for 
every s ^ d, if p is a random state on distributed according to pd,s, o,nd if X is 
a standard Gaussian vector in Hq, then 



E 



Applying the lemma for K 



I 

d ^ 
AWTo 



1 



E\\X\\k. 



d^/s 

we obtain that 



E 



P- 



AVVTo 



So{d) 



This suggests that the threshold for the set AWTq occurs at So{d), since a state p 
is APPT when \\p - l/d\\_AvvTo < 1 and non-APPT when ||p - l/d\\AvvTo > 1- 

The proof of Theorem 3.3 is almost identical to that of Section 4 in [4], and here 
we sketch its proof for completeness. We refer the readers to [4] for more details, in 
particular Appendix E for the Levy's Lemma and concentration of measure theory. 

Proof of Theorem 3. 3. Let Shs be the Hilbert-Schmidt sphere in the space of d x s 
matrices (it can be identified with the real sphere 5'^'^'^"^) and / : Shs — )■ R be the 
function defined by 

/(M) = MM* - ^ 



AWT, 
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Formula (11) asserts that E/ ^ 




. For every r > 0, denote by Q = Q{r) the 



subset 



n = {MeSHs : ||M|U^r}. 
Inequahty (4) and the proof of Lemma 4.2 in [4] imply that the Lipschitz constant of 
/|q is bounded by 2rd. Note that Q{r) = Shs since ||M||oo ^ ||^||2 = 1- Then, the 
global Lipschitz constant of / is bounded by 2d, and hence the median of / (denoted 
Mf) differs from its mean, E/, by at most C2d/\/2ds = C \fdfs (see Appendix E 
in [4] ) . It follows that the median of / is also of order ^J sq/ s. 

By a net argument similar to that in [3], one has ¥{Shs \ ^) ^ exp(— cs) if 
r = "i/y/d. A local version of Levy's lemma (see Appendix E in [4]) implies that for 



whenever Mf ^ 2/3 (or, equivalently, s > sq). This ends the proof of Theorem 



We consider the product systems = C^^ ® C^^. Recall that a random state p 
on = C^^ (g) C^^ distributed according to fid,s has the same distribution as MM* , 
where M is a. d x s matrix uniformly distributed on the Hilbert-Schmidt sphere of 
the Hilbert space of ci x s complex matrices. A more convenient, but equivalent, 
way is to link the measure fid,s with a normalized Wishart matrix. More precisely, 
let W = GG\ where G e M^ xs(C) is a Ginibre matrix, i.e. a matrix with i.i.d. 
standard complex Gaussian entries. The random matrix W is called a Wishart 
matrix of parameters (rf, s). A random induced state p distributed according to pd,s 
is then given by the formula p = W/tiW [25, 44]. 

The following theorem is the main result of this section, providing asymptotic 
values of the threshold for the convex body AWT. This theorem covers the case 
when p = mm{di,d2) — ?■ oo. The simpler case of p bounded is treated in Theorem 
4.2 as it gives sharper bounds. 




3.3. 



□ 



4. Estimates on the threshold for AWT 
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Theorem 4.1. Let p be a random induced state distributed according to the proba- 
bility measure fid,s- 

(i) For all e > 0, almost surely, when d ^ oo and s > (4 + e)p^d, the quantum 
state p is APPT; 

(ii) When 1 <^p^ <^ d and s < (4 — e)p^d, p is not APPT almost surely; 

(iii) When ~ rd for a constant r G (0,1], there exists a constant Cr (see 
formula (13) ) such that whenever and s < 4(Cr — e)p'^d, p is not APPT 
almost surely. 

Proof. We start with (i). For given eigenvalues {Ai, ■ ■ ■ , A^}, we introduce the fol- 
lowing p X p matrix: 

T = T(A; a+, (t_) = e(A; (t+, a_) - 2d-\, 

where Ip denotes the p x p identity matrix. Recall that p = min((ii, ^2). From 
formula (6) and Theorem 2.7, the matrix T has small entries: |Tjj| ^ (4 + e)/\fds. 
A necessary condition for the matrix 6 = 2d~^lp + T to be semi definite positive is 
that T should have operator norm smaller than 2/d. It is a well known fact in matrix 
analysis (see [21]) that 

||T||op ^ p||T||i_,oo = P max ^ ^j}^ . 

The conclusion in the statement follows by asking that (4 + e)p/\fds ^ 2/d. 

We move now to the proofs of (ii) and (iii). We shall proceed by exhibiting a vector 
X G R^, such that, X"^A(A; cr+, a^)x < for some pair of linear orderings. This does 
indeed suffice to show that the matrix O is not semidefinite positive. Indeed, we take 
the column vector a; = (l,l,...,l)'^G R^. Any pair of linear orderings is compatible 
with such a vector, and one has 

p p+ p- 

(12) x'A(A; a+, cr„)x = ^ Aij = ^ Xd+i-i - ^ Ai, 

i,j=l i=l i=l 

where p+ = p{p + l)/2 and p_ = p{p — l)/2. 

We shall now consider the two regimes in the statement, starting with 1 <^ <^ c?. 
The main idea here is to note that, for all z, Aj = l/rf + Aj/A/ds, where Aj are the 
eigenvalues of the matrix introduced in Section 2.2. By Theorem 2.7, for all e > 
and for d large enough, all the "large" eigenvalues Ai, . . . , Xp_ appearing in equation 
(12) are bigger than l/d+ {2 — e)/\/ds; in the same vein, all the "small" eigenvalues 
Xd+i-p^ , . . . , Arf are smaller than 1/d — {2 — e)/ \fds. We obtain 

X A(A; (T+, a^)x ^ P+ 3 j=- - 3 



d Jds j \d Jds ) d Jds 
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which is negative as long as s < (4 — e)p^d. This concludes the proof for the first 
regime 1 d. 

We now move on the the second regime, where ~ rd, for a fixed constant 
r G (0, 1]. Writing equation (12) in terms of the eigenvalues of Zd, we obtain 



p 1 
x'A{X; 0-+, a-)x = - + -= 
" Vds 



Using the fact that the asymptotic spectrum {A.j} of is semicircular, we obtain 
the following bounds for d large enough: 



^ Ad+i_i < / xw(x)dx, and 

i=l >^-2 ^ i=l 



where 

u'(a;) = — V4 — 

is the density of the standard semicircular distribution and ^ [0,2] is defined 
implicitly by 

f2 

w(x)dx = r/2. 

Indeed, it is a classical result in random matrix theory that the conclusions of 
Theorem 2.7 and of Corollary 2.5 for the matrix model Z^ imply that its repartition 
function converges almost surely uniformly towards the repartition function of the 
semi-circle distribution. This implies the previous claim. 

By the previous bounds, we obtain that 



2 

xw{x)dx. 



Using p ~ -y/rd^/^, the above expression is seen to be negative as soon as 



, xw(x)dx , 
s<^fd\'-^ =A{C^-e)p'd, 



2 



where we put 

(13) a 



xw{x)dx 



Cr/2 
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One can easily show that the map r h-> C,- is increasing and that 

Cq = hm = 1 and C*i = I / xw(x)dx ] — , 
r->o \Jq ) Qtt^ 

and the proof is complete. □ 

When p = min((ii, ^2) is fixed and s/d ^ c for a constant c > as — 00, one 
can obtain the following sharp estimate on the threshold for AWT. 

Theorem 4.2. Let p he a random induced state distributed according to the measure 
Pd,s- Almost surely, when d ^ 00 and s ~ cd, one has: 

(i) p e AWT, ifc>{p+ Vp2^)2. 

(ii) p i AWT, ifc<{p+ ^/p^3I)2_ 

To prove this result, we need the following well-known result in random matrix 
theory. This result describes the behavior of the spectrum of a Wishart matrix of 
parameters {d,s) with s/d ^ c and d ^ 00 [6]. 

Proposition 4.3. Let Ai ^ A2 ^ • ■ ■ ^ A,^ ^ 6e the eigenvalues of a Wishart 
matrix of parameters {d, s) . Then, in the asymptotic regime d — )■ 00 and s ^ cd for 
a constant c > 0, one has: 

(i) Almost surely, when d 00, the empirical eigenvalue distribution 

d 

1=1 

converges weakly to the Marchenko-Pastur ( or free Poisson ) distribution tt^ 
given by the formula 



AlC — (x — 1 — c)^ 
TTc = max(l - c, 0)(5o H l[(^„i)2 (^+i)2](x) dx\ 

(ii) For every fixed integer k, almost surely, as d 00 

'0 ifc^l, 



Ad, Arf_i, . . . , Xd^k+i o-c 
and 



Ai,A2,...,Afe->6e = (yc+l)l 

Proof of Theorem 4-2. Recall that a spectrum (Ai ^ ■ ■ ■ ^ A^) corresponds to states 
in AWT if and only if the matrix 0(A; (t+, cr_) = A(A; 0"+, (t_) + A(A; 0"+, cr„)^ is 
positive for all pairs (cr+, cr_). Such a criterion for APPT is invariant under scaling of 
the matrix p. Thus, it is equivalent to consider the non-normalized Wishart matrix 
W. 
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In our case, we have that, for all i < j, a^{i,j) ^ di{di + l)/2 and a^{i,j) ^ 
di{di — l)/2, which are bounded quantities. Hence, it follows from Proposition 4.3 
that, asymptotically, matrices 6(A;(T+,(T_) are all equal to 



A. 



2ac ac-bc ■ ■ ■ ttc- &c\ 
flc - be 2ac ■■ ■ ttc- b, 

\ac - be ac - be ■ ■ ■ 2ae / 



+ be)Id, + {ae-be)l 



di 



The eigenvalues of the matrix above are (oc + be) + c?i(ac — be) with multiplicity one 
and (oc + be) with multiplicity di — 1. 

Note that if c ^ 1, then ae = and be > 1 which implies that {ac+be)+di{ae — be) < 
0. Equivalently, the matrix is not positive and (Ai, ■ ■ ■ , A^) does not correspond 
to states in AWT. On the other hand, if c > 1, Ac is positive if and only if 
+ be) + di{ae — 6c) ^ which can be shown to be equivalent to the condition in 
the statement. □ 

Combing Theorems 4.1 and 4.2, we get the following theorem; although this state- 
ment is strictly weaker than the results above, it captures the behavior of the thresh- 
old in an unique statement. 

Theorem 4.4. There are effectively computable absolute constants c, C > 0, such 
that, if p is a bipartite random induced state on = C^^ ® 0"^^, obtained by partial 
tracing a random pure state on ® C^, then for d = did2 large enough, 

(i) The random density matrix p is not AWT with very large probability when 
s ^ cp'^d; 

(ii) The random density matrix p is AWT with very large probability when s ^ 
Cp'^d. 

The above theorem asserts that the thresholds for AWT is indeed (approxi- 
mately) 4p'^d. Together with Theorem 3.3, one can obtain the estimate for the mean 
width of AWT°. 

Corollary 4.5. Let AWT be the set of states with APPT on the bipartite system 
Qd _ Qdi g>, Qd2 ^ j'liQ^ ifiQ threshold function Sq for AWT satisfies 

So = so{di, d2) ~ w{AWT°)'' ~ p'd. 

In particular, w(AWT°) ~ py/d. 
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